Modified KdV hierarchy : Lax pair representation and bi-Hamiltonian structure 
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We consider equations in the modified KdV (mKdV) hierarchy and make use of the Miura trans- 
formation to construct expressions for their Lax pair. We derive a Lagrangian-based approach to 
study the bi-Hamiltonian structure of the mKdV equations. We also show that the complex modi- 
fied KdV (cmKdV) equation follows from the action principle to have a Lagrangian representation. 
This representation not only provides a basis to write the cmKdV equation in the canonical form 
endowed with an appropriate Poisson structure but also help us construct a semianalytical solution 
of it. The solution obtained by us may serve as a useful guide for purely numerical routines which 
are currently being used to solve the cmKdV eqution. 
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1. Introduction 

Nearly forty years ago Lax [1] showed that the Ko- 
rteweg de Vries (KdV) initial value problem for u — 
u{x,t) given by 



with 



u{x,0) = Vix) 



(1) 



(2) 



is but one of the infinite family of equations that leave the 
eigenvalue of the Schrodinger equation with the potential 
V{x) invariant in time. The subscripts of u in (1) denote 
differentiation with respect to the associated independent 
variables. The family of equations discovered by Lax 
often goes by the name KdV hierarchy and is generated 
by making use of the recursion operator [2] 



4u 



OX 



in the differential relation 



0, 1, 2, 3, 



(3) 



(4) 



The KdV equation (1) is recognized as the solvability 
condition for the system 



and 



with 



dtip = Atp ,dt 



L = -dt 



d_ 
dt 



(5a) 



(55) 



(6a) 



the so-called Schrodinger operator, 
order linear operator written as 



Here A is a third- 



(66) 



The existence of the solution i(j 
constant A is equivalent to 



V'(A, X, t) for every 



dtL = AL- LA= [A, L] 



(7) 



The result in (7) is called the Lax equation and the op- 
erators L and A are called Lax pair [1]. The Lax pair 
representation holds good for all equations in the KdV 
hierarchy. In the context of Lax's method it is often said 
that L defines the original spectral problem while A rep- 
resents an auxiliary spectral problem. As one goes along 
the hierarchy, L remains unchanged but the differential 
operator associated with the auxiliary spectral problems 
changes according to 



A^^ 



0, 1, 2, 3,, 



(8) 
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The operator Aq = dx and Ai stands for A in {6b). It 
is not always easy to obtain results for other A^s which 
generate higher-order KdV equations. The coefficient aj 
depends on the solution u and derivatives Un{— f;pr)- 
From (8) it is clear that as j varies, the dimension of aj 
changes. Thus aj should be chosen as a linear combina- 
tion of power and products of u and u„ 's such that the 
terms in the curly bracket have the right dimension of 
g2n+i^ The constructed expression for An will then gen- 
erate the KdV hierarchy when used in the Lax equation 
[3]. On the other hand, one can postulate that for an 
evolution equation of the form Ut — K[u] the terms in 
the Frechet derivative of K[u] contribute additively with 
unequal weights to form the operator An such that L 
and An via (7) reproduces the equations in the hierarchy 
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[4]. Of course, there should not be any inconsistency in 
determining the values of the weight factors. 

Zakharov and Faddeev [5] developed the Hamiltonian 
approach to integrability of nonlinear evolution equations 
in one spatial and one temporal (1+1) dimensions and, 
in particular, Gardner [6] interpreted the KdV equation 
as a completely integrable Hamiltonian system with dx 
as the relevant Hamiltonian operator. A significant de- 
velopment in the Hamiltonian theory is due to Magri [7] , 
who realized that integrable Hamiltonian systems have 
an additional structure. They are bi-Hamiltonian i.e. 
they are Hamiltonian with respect to two different com- 
patible Hamiltonian operators dx and (i9^ — 4,udx — 2ux) 
such that 

n = l, 2, 3... . (9) 

Here H,i = j Tindx with H,, , the conserved densities for 
the equations in the KdV hierarchy. These conserved 
densities generate flows which commute with the KdV 
flow and as such give rise to an appropriate hierarchjf. 
Traditionally, the expression for 7i„ is constructed using 
a mathematical formulation that does not make explicit 
reference to the Lagrangians of the equations in the hi- 
erarchy. However, a Lagrangian-based approach can be 
used to identify 7Y„ as the Hamiltonian density of the nth 
hierarchical equation [8]. 

The nonlinear transformation of Miura or the so-called 
Miura transformation [9] 

u = Vx+v'^, v = v{x,t) (10) 

converts the KdV equation into a modified KdV (mKdV) 
equation 

Vt = V3x - Qv'^Vx . (11) 

This equation differs from the KdV equation only be- 
cause of its cubic nonlinearity. It has many applicative 
relevance. For example, mKdV equation has been used 
to describe acoustic waves in anharmonic lattices and 
Alfven waves in collisionless plasma. It is of interest to 
note that the recursion operator A for the mKdV equa- 
tion [10] 

= dl - iv^ - 4vxd-\v (12) 

can be identified from (4) and (10). The equation of the 
mKdV hierarchy can be generated by using the relation 

t;^ = A>, , n = 0, 1, 2, 3,... . (13) 

It is straightforward to obtain the equations in the mKdV 
hierarchy from those in the KdV hierarchy by the use of 
Miura transformation. However, it is a nontrivial prob- 
lem to derive the Lax representation and construct the 



bi-Hamiltonian structure of the equations in mKdV hi- 
erarchy starting from corresponding results for the KdV 
equation [11]. In this work we shall deal with these prob- 
lems. To derive the Magri structure we shall make use of 
a Lagrangian-based approach. 

In addition to the above another system of our interest 
is the complex modified KdV (CMKdV) equation given 
by 

Vt = V3x - Q\vfVx , (14) 

This equation follows from the third-order nonlinear 
Schrodinger equation via an appropriate variable trans- 
formation [12]. We shall provide a variational formula- 
tion of (14) which, on the one hand, allows us to study 
its canonical structure and, on the other hand, serves 
as a useful basis to construct an approximate analytical 
solution in term of a trial function. In this context we 
note that the numerical routine for solving such equa- 
tions is quite complicated [13] and requires the use of 
Crank Nicolson method for time integration and quintic 
B-spline function for space integration. We believe the 
solution presented by us may serve as an initial guide for 
the more ambitious programmes. 

In §2 we introduce the equations in the mKdV hi- 
erarchy and derive their Lax pair representation. We 
find that the system of equations follows from the ac- 
tion principle and as such can be obtained from appropri- 
ate Lagrangian densities via the so-called Euler-Lagrange 
equations. The corresponding Hamiltonian densities con- 
stitute the involutive conserved densities of the mKdV 
equation. We then study the bi-Hamiltonian structure 
of the mKdV equations. In §3 we convert (14) to a vari- 
ational problem and thus obtain a Lagrangian represen- 
tation for the equation. As a useful application of the 
Lagrangian density so derived we work out the canoni- 
cal form [5] of the cmKdV equation and also construct 
a solution of it by means of sech trial functions and a 
Ritz optimization procedure [14]. Finally, in §4 we try to 
summarize our outlook on the present work. 

2. mKdV hierarchy 

The equations of the mKdV hierarchy follow from 
(13) for n = 0, 1, 2, 3... . We shall construct Lax pair 
representations of these equations by taking recourse 
to the use of (10) in (6a) and (8). For these equations 
we shall use a Lagrangian-based method to obtain the 
conserved densities which are in involution and generate 
the so-called mKdV fiow. We shall then try to realize the 
bi-Hamiltonian structure by an appropriate modification 
of (9) by the use of Miura transform. 

(a) Lax pair representation 

From (6a), (66) and (10) we write 

L = -dl + v^+ Vx (15a) 
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and 

A = 4dl- 3{v^ + v^)d^ - 39^.(^2 + v^) . (156) 
Using (15) in (7) we get 

{d^ + 2v) {vt - + Gv'^v^) = . (16) 

As with (11), (16) gives the mKdV equation. In view of 
this we shall denote the Lax pair in (15) by and A™ 
just to indicate that these refer to the mKdV equation. 
We shall follow this convention for all operators and func- 
tions related to the mKdV equation. Consistently with 
the notation of (8) A™{= A) stands for AJ'-. In close anal- 
ogy with the case of higher KdV equations the original 
spectral problem for the mKdV equations characterized 
by the operator remains unchanged as we go up the 
hierarchy but the differential operator ^™'s change with 
n. From (10) and the results given in refs 3 and 4 one 
can calculate the expressions for A™ , n = 2, 3, 4, ... . In 
the following we present some of our results. 

A^ = 16dl + {25vl + SOv'^v^ + 10vv2a= + I5v^ + S^^sc.)^^- 

20dl.{v^ + v^) , (17a) 
= 64dl - UO{v* + 3vl + 2v\^ + 2vv2^ + t;3x)5^+ 
112(t;2 + v^)dl^ 140^1 («'' + Zvl + 2v'^v^ + 2vv2^+Vi^)+ 
\Vldl.{v^ t;^) (70u'' 2V^v\^ -|- 1050t;\^ -|- 210t;^-F 
\^<^v^vix + M<^vVxV2x + 721u^^ -h l^v^vzx)dx^ 

d^.{70v^ + 210v^v^ + W50v\l + 210vl) 

+dx.{U0v^V2x + 840vVxV2x + 721vl^ + 70v'^V3x + 

798t;xt;3x + 1821;^^" + 91v5x) (176) 

and 

AX = 256dl-255vsx-'i-79^V7xdx-510vvrx-1152vxdl- 
1152u^a^-5628i;6x5^+2100t;xi'6x-3588i;u6x9x+5670v\6a;- 

4032t;2^5^ - 8064vVxd^ - 10248z;5,.a3 - 7224i;2^i;5r.- 



11256t;t;5x9^-15594t;^V5a;5^+18312t;t;^t;5^-|-5934u2v5^a^-|- 
11340v3v5a; - 8736v3xd^ - Il'il2vv2xdl - lMhQvldl+ 
AQ32v\^dl + 2QlQv'^dl - 11760u4^a^ - 11760v3^V4^- 
2M9&vvixdl - 392Mv2xVixdx + 19152t;u2a:t^4x- 
A3%mvxVAxdl+l2mQv'^Vixdl+&8Q7QvlvAx+^^'^^^vvxVAxdx+ 
70224t;^u^t;4a: + 11868t;^t;4a:9x - 210t;%x - 23520ut;3^a^- 

m2AQvxV2xdt+^Q32Qv\2xdt+'2'Q&'^'^vvldt+2QQA0v\xdt- 
25758v|^a^ + 12180i;u|^ - 66864u^i;3^5^ + IbUQv^vzxdl- 

873Q0v2xV3xdl + nQ2Q8vxV2xVzx + Qi972QvxV2xV3xdx+ 

13Q2Q(iv%2xVsx + 7^&QQvVxV5xdl + 2^2QQv^vsxdl+ 

81660u^U3x9x + 64596ut;^t;3a; + 93336u^Ua;U3x9a;- 

IhQmv^VxVzx - QZmv'^vrixdx - ^2Qv^V5x - 50568w2a:i9;^+ 

7392QvVxV2xdl+28hmvldl+Q888Qv\ldl-hQAQv'^Vxdl- 

1680t;^5^ + 19026t;|^ + bMmvvl^dl + llAAmvxvl^dx+ 

m'ih2vVxvl^+Q7272v\l^dx-lhl2Qv^vl^+llMAQvlv2xdl+ 

l&l28Qv\xV2xdl - 7bQQv'^V2xdl + 208488i)u^t;2x9x+ 

57960t;^U2x - 66780i;^t;^i;2x - &(iA8Qv^VxV2xdx- 

27720i;^?;xi;2x - l2mQv^V2xdx + 1260i;^i;2x+ 

85680vt;^a^ - 30240u^t;^a^ - l^l2Qv^Vxdl + 28518t;^ax- 

27720t;t;^ - 57960t;^t;^ax - 378001;^!;^ - AAlQQv'^vldx+ 

7bmv^vl + 2b2Qv^v^dx + 2b2Qv^Vx + Q^Qv^dx ■ (17c) 
Using (15a) and (17a) in (7) we get 

{d^+2v){vt-V5^+A0vv^V2x+^0v\3^+10vl-30v^v^} = . 

(18) 
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The expression inside the curly bracket represents the 
equation obtained from (13) with n — 2. Results similar 
to that in (18) hold good for any pair like [A™, L™]. 
This observation serves as a useful check on our results 
for A™ with arbitrary values of n. 



(256) 



= \wtw^ - ]^wl, -wl- bwlwl, , (24c) 



(b) Bi-Hamiltonian structure 

Here we shall demonstrate that the bi-Hamiltonian 
structure of (11) and all higher-order equations obtained 
from (13) with n = 2, 3, 4.... We note that a single 
evolution equation hj = P[u], u e R is never the Euler- 
Lagrangc equation of a variational problem [10]. One 
common trick to put a single evolution equation into a 
variational form is to replace u by a potential function 



V = —Wx, w = w{x,t) 



(19) 



(25c) 



11 7 



35 



70 4 2 5 



(24rf) 



The function w is often called the Casimir potential. Our 
expressions for the Lagrangian densities will be written 
in terms of w and its appropriate derivatives. Hamil- 
tonian densities obtained by use of Legendre map can, 
however, be expressed in terms of field variable v{x,t) 
and its derivatives. 

The linear equation obtained from (13) with n = 
reads 



Vt = Vx 



(20) 



1 7 

Ti-T = -VV6x - -^V^Vix - l^V^VxVSx 



21 



2 2 
' 'V2X- 



35 2 35 c TO 4 2^8 / T\ 

YVVxV2x + y V V2x - Y'" 2^ ^ ' 



and 



= ^WtWx - ^WxWQx + ^wlwrx + 27wlw2xW6x + 



From (19) and (20) 

Wxt = Wxx = P[wx] (say) 



(21) 



^r, 2 105 2 ^1 >s 69 2 2 

57WxW3xW5x + -^WxW2xW5x " ^i-WxW^x + y ^^^4^ + 



Equivalently, 



Wt = Wx 



(22) 



91 

189wxW2xW3xW4x - 168wlw2xWix + y w^xM'L - 126w^wL 



In writing (22) we have used the boundary condition limit 
w{x,t) = as a; ^ ±00. The self adjointness of Plwx] 
ensures the existance of a Lagrangian for (21) and (22). 
In this case, the Lagrangian density can be constructed 
using the homotopy formula [10] 



m= [\p[K]dX- (23) 
Jo 



From (23) we get 



105 315 
-518wlwlxW3x+^wlw3x-'i-3Swlwix+-^wlwlx-7wl° 

(24e) 



1 9 
-—1 

2'""" 2 



105 

Hf = ^vv8x-7;V^vex-'^7v'^VxV5x-57v'^V2xV4x — ^vvlvix 



69 91 

+21v^V4,x —v'^vl^-189vVxV2xV3x + '^Q8v%xV3x ^'"'"2x 



^WtWx - -Wx 



(24a) 



The subscript zero is self explanatory. The Hamiltonian 
density obtained from (24a) is given by 



I-Lq 



1 2 

2^^ 



(25a) 



The Lagrangian and Hamiltonian densities for the mKdV 
(n = 1) and higher-order equations obtained from (13) 
for n = 2, 3 and 4 are given by 



£5" 



1 



1 



zWtWx 



:WxW3x 



1 4 

2"^x, 



(246) 



105 

+126v\lx + ^18v^vlv2x - -i^v''v2x+ 



133^2^4 _ 315^6^2 ^ 7^10 _ 



(25e) 



Results of WJ^'s for still higher values of n can be ob- 
tained in a similar manner. As a useful check on our 
expressions, one can verify that these results are in ex- 
act agreement with those obtained by the application of 
Miura transformation on the well known conserved den- 
sities of the KdV equations. 
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The bi-Hamiltonian structure of equations in the 
mKdV hierarchy can easily be verified by using our 
Hamiltonian functionals in 



where 



£={dl-Av^d,-4v,d-\vd,) 



(26) 



(27) 



The first Hamiltonian operator dx in (26) is the same as 
that in (9) while the second has been obtained from [10] 



(28) 



The operators dx and £ are skew symmetric and satisfy 
the Jacobi-identity. Thus they constitute two compatible 
Hamiltonian operators such that all equations obtained 
from (13) are integrable in the Liouville's sense [7]. 

3. cmKdV equation 

The complex mKdV equation in (14) can be restated 
as a variational problem given by 



dt \ dvt J 6v 



d f dC\ SC 



dt \ dvl ) Sv* 



with the variational derivative 

3 



= 



5tp 



n>0 



dipn 



Here 



(316) 



(32) 



(33) 



It is easy to verify that (30) and (31&) give the cmKdV 
equation while (30) and (31a) give the corresponding 
complex conjugate equation. The Hamiltonian corre- 
sponding to the Lagrangian density (30) is given by 



I' 



H = / ndx 



(34) 



with the Hamiltonian density 
1 3 

In order to show that (14) is a Hamiltonian system we 
will have to write it and its complex conjugate in two 
different forms 



and 



vt = {v*{x) , H{y)} 



v; = -Mx) , H{y)} . 



(36) 



(37) 



We have already found an expression for the Hamilto- 
nian. Thus our task is to look for fundamental Poisson 
bracket relation for the field variables that reduce (36) 
to the cmKdV equation and (37) to the complex conju- 
gate one. One can easily check that the required Poisson 
bracket relations are given by 



{v{x) , v{y)} = {v*{x) , v*{y)} = 5{x - y) 



(38) 



The relations (36) and (37) can be written in the sym- 
plectic form 



5 j J ^{v, V*, Vx, vl, V3x, v^x, Vt, Vt, X, t) dx dt = 

(29) 

with the Lagrangian density written as 

11 3 

= 2 i^*'"t - ^0-2 (^*^3^ ~ vv3x)+-^vv* {v*Vx - vv*) . 

(30) 

The Euler-Lagrange equations corresponding to (29) are 

(31a) 
and 



with J = 



Vt 

1 
-1 



or] 



(39) 



a skew-adjoint matrix as the 



Hamiltonian operator. 

Equation (14) arises in a number of applicative con- 
texts including the nonlinear evolution of plasma waves 
[15]. To our knowledge there is no well-defined spec- 
tral problem that can easily be used to solve the cmKdV 
equation in terms of known transcendental functions. 
But a number of works has been envisaged to obtain the 
solitorjf- waves and/or soliton solutions of this equation. 
See, for example, [13] and references therein. Here we are 
interested to provide an accurate approximation solution 
of (14) by supplementing the Lagrangian density in (30) 
with sech trial functions and a Ritz optimization proce- 
dure. We have chosen to work with the trial function 
written as 



v{x, t) = a{t)sech [{x - y{t))/'w{t)] x 



[i{q(t)+r{t)(x-y(t)) + ^{x-v(t)f)] 



(40) 



Here the parameters a, y and w are related to the three 
lowest-order moments of the v envelope and represent, re- 
spectively, its amplitude, central position and width. The 
other parameters q, r and b stand for the phase, velocity 
(center of mass) and frequency chirp. Understandably, 
these parameters will all vary with time t. Using (40) in 
(30) we get 



(41) 



6 



where 



we obtain 



w 



w 



+ 



(46a) 



w 



a ry seen ( \ — a q seen 



w 



-4au>r — ISwa r — it ab rw w ab — 4au'g+ 

w 6 



a bwsech ' 



2 \ w 



(42a) 



Aawry + —abtuw = 
6 



(466) 



£(2) = ^,ec/.^(^han/i^^"-^ 



3a^r 



sech 



w 



w 



x-y 
w 



a^r^ seeh^ 



2 I x-y 
w 



x-y 
w 



a?b'^r seeh^ 



2 / x-y 
w 



and 



£(3) =-3aVsec/i^ 



4 I x-y 
w 



(426) 



(42c) 



Here the dots stand for derivative with respect to t. The 
subscript s on £ merely indicates that we have inserted 
the sech ansatz for v{x, t) into the Lagrangian density. 
In terms of (41) the variational principle (29) leads to 

5 j <L> dt = , (43) 

with the averaged effective Lagrangian 

<L>= Lsdx . (44) 

J 00 

The result for < L > is given by 



r) o A ''0 

< L >= —2wa r — Awa r — —a b rw 



2ofr 



w 



—w^a^b — 2a^wq + 2a^wry + —a^bww . (45) 

The reduced variational principle expressed by (43) re- 
sults in a set of coupled ordinary differential equations for 
the parameters of our trial function. From the vanishing 
condition of the variationals 

S<L> S<L> S<L> S<L> S<L> 



6q 



5a 



and 



6y 



5<L> 
Sb 



6w 



6r 



-2d^r' - Aa^r - '—af'h^r + 



2df'r TT" 



(46c) 



wa^b - 2arq+ 



7^2 ^ /^2 \ 

2a'ry + -^a%w - - I -^a%w \ = , (46d) 



2a^ 



-6«;a^r^ - 4«;a'' - —a^Vw h 2a^wy = (46e) 

2 w 



and 



-T,^a%rw + —a^w + ^ f —w^a" 1=0. (46/) 
12 dt \12 ' ^ •'^ 



Equations in (46) can be used to write 

a^w — constant = Eq , 



r = constant , 



3a6r 



w 



da 
'dt 

^=3r^ + 2a^ + 4b' + ^, 
dt 



dw 



6br 



and 



db _ 24r a^ 24r 1 

dt TT^ W TT^ 



(47a) 



(476) 



(47c) 



(47d) 



(47e) 



(47/) 



Equation (47a) expresses the variational version of the 
energy conservation law [16] while (476) states that the 
center of mass of the solution of (14) moves with a con- 
stant velocity. For a given values of r, the set of coupled 
ordinary differential equations (47c)- (47/) can easily be 
solved numerically. Note that knowledge of a(t), y{t) and 
w{t) can be used to study the \v{x, t)\ as functions of x 
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FIG. 1: |v(a;,t)| as a function of t for three different 
values of x. Here r = 0.1. 
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FIG. 2: Same as that in FIG. 1 but with r = 0.001 



realize the bi-Hamiltonian structure. 

In close analogy with the mKdV equation, the cmKdV 
equation in (14) also follows from Hamilton's variational 
principle provided the action functional is made to 
vanish for simultaneous variations of both v and v* . In 
this case the Lagrangian density is a function of v, v* and 
their appropriate time and space derivatives. We could 
use the Hamiltonian corresponding to this Lagrangian 
density to write the cmKdV equation in the canonical 
form [5] with an appropriate Poisson structure. As 
an added realizm we demonstrate that the Lagrangian 
density constitutes a basis to derive a semianalytical 
solution of (14). We achieve this by taking recourse 
to the use of sech trial functions to define a reduced 
variational problem which in conjunction with the Ritz 
optimization procedure could yield an uncomplicated 
solution of the cmKdV equation. There exist some 
sophisticated numerical routines [13,15] to solve the 
equation. However, we feel that the variational approach 
sought by us will serve as a complementary tool towards 
understanding the properties of solitory wave- and/or 
soliton-solutions of the complex modified KdV equation. 
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and t. We worked with the initial conditions a(0) = 1, 
&(0) = 0, 1/(0) = 0, w{0) = 1 and solved these equations 
using the fourth-order Runge-Kutta method [17]. First 
we take r = 0.1 and plot in FIG.l, \v{x, t)\ as a function 
of t for three different values of x, namely, a; = (solid 
curve), a; = 10 (dotted curve) and x = 20 (dashed curve). 
From these curves it is clear that as x increases \v{x, t)\ 
decreases rapidly. This implies that for our chosen value 
of the velocity we have decaying solitory wave solution. 
We have verified that for still higher values of r the solu- 
tions decay more rapidly. In FIG. 2 we present a similar 
plot of \v{x, t)\ for r = 0.001. Interestingly, \v{x, t)\ re- 
mains unchanged as x increases. Thus one can infer that 
the solutions of (14) for small values of r behaves like 
solitons. 

4. Conclusion 

The nonlinear transformation of Miura or the so-called 
Miura transformation is an aid to obtain the modified 
KdV (mKdV) equation from the KdV equation. We find 
that this transformation also provides an effective way 
to construct expressions for Lax pairs of all equations in 
the mKdV hierarchy. As with the KdV equations the bi- 
hamiltonian structure of the mKdV equations are tradi- 
tionally studied using involutive set of conserved Hamil- 
tonian densities without explicit reference to their La- 
grangians. We derive a Lagrangian-based approach to 
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